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Abstract 

Background: This study is aimed at the analysis of genetic and physiological effects of myostatin on economically 
relevant meat quality traits in a genetic background of high muscularity. For this purpose, we generated G 3 
populations of reciprocal crosses between the two hypermuscular mouse lines BMMI866, which carries a myostatin 
mutation and is lean, and BMMI806, which has high intramuscular and body fat content. To assess the relationship 
between muscle mass, body composition and muscle quality traits, we also analysed intramuscular fat content 
(IMF), water holding capacity (WHC), and additional physiological parameters in M. quadriceps and M. longissimus in 
308 G 3 -animals. 

Results: We found that individuals with larger muscles have significantly lower total body fat (r = -0.28) and IMF 
(r = -0.64), and in females, a lower WHC (r = —0.35). In males, higher muscle mass was also significantly correlated 
with higher glycogen contents (r = 0.2) and lower carcass pH-values 24 hours after dissection (r = —0.1 9). Linkage 
analyses confirmed the influence of the myostatin mutation on higher lean mass (1.35 g), reduced body fat content 
(-1.15%), and lower IMF in M. longissimus (-0.13%) and M. quadriceps (-0.07%). No effect was found for WHC. A 
large proportion of variation of intramuscular fat content of the M. longissimus at the myostatin locus could be 
explained by sex (23%) and direction-of-cross effects (26%). The effects were higher in males (+0.41%). An 
additional locus with negative over-dominance effects on total fat mass (-0.55 g) was identified on chromosome 
16 at 94 Mb (86-94 Mb) which concurs with fat related QTL in syntenic regions on SSC13 in pigs and BTA1 in 
cattle. 

Conclusion: The data shows QTL effects on mouse muscle that are similar to those previously observed in 
livestock, supporting the mouse model. New information from the mouse model helps to describe variation in 
meat quantity and quality, and thus contribute to research in livestock. 



Background 

In livestock production, there is a high interest in control- 
ling meat quantity and quality; knowledge about genes 
affecting muscle size and other meat properties can help 
breeders to select animals according to desired traits. 
Myostatin (Mstn, or growth and differentiation factor 
8 - Gdfi) was first identified in knock-out mice as being a 
gene responsible for regulation of muscle growth [1]. 
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However, the hypermuscular phenotype that originates 
from a hyperplasia and hypertrophy of muscle fibres was 
known long before its molecular-genetic background was 
elucidated [2]. In sheep, pigs and cattle, several mutations 
of myostatin or its promoter region have been identified 
as affecting muscle size [3-6]. Furthermore, myostatin also 
influences glucose metabolism and fat accumulation, as 
shown in knock-out mice that had smaller adipocytes [7] 
and did not develop obesity on a high-fat diet. These 
results suggested an altered metabolism for the utilization 
of lipids and glucose as energy fuel that prevented insulin 
resistance in mice [8,9]. A decreasing effect on intramus- 
cular fat (IMF) content and carcass fat proportion was also 
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described for the "double-muscled" phenotypes in cattle 
[10,11], which was later found to be caused by a myostatin 
mutation as well. However, different myostatin mutations 
can cause different degrees of hypermuscularity. The most 
extreme form of this phenotype in cattle is seen in the 
Belgian Blue, while different mutations with less extreme 
phenotypes were identified in other breeds [12]. 

Due to its large implication, not only on skeletal muscles, 
a lot of research has been undertaken on the effects of 
myostatin and its possible applications in health care and 
livestock production. The myostatin-gene, for instance, 
appears to be a promising candidate for the treatment of 
muscle dystrophy diseases [13,14], cardiac tissue regener- 
ation or metabolic syndrome [15,16]. In livestock, efforts 
to manipulate the gene's expression through immunization 
have been made for the purposes of higher lean mass 
production [17]. In accordance with this, modifiers of 
myostatin, such as follistatin (Fst), which inhibits myosta- 
tin, are currendy under review as therapeutics [18,19]. 

In the present study, we generated G3-populations of 
reciprocal crosses between two hypermuscular Berlin 
Muscle Mouse Inbred (BMMI) lines to examine the gen- 
etic characteristics of myostatin and to find additional 
genes that could influence muscle growth and compos- 
ition. The BMMI lines were hypermuscular as a result of 
long term selection that mirrors the selection process for 
high meat yield in livestock. The parental BMMI866 line 
originates from a population with the "compact" pheno- 
type and therefore carries the known Mstn Cmpt ' dllAbc mu- 
tation [20,21]. Although the BMMI806 line originates 
from the same founder population, it does not carry this 
mutation; BMMI806 animals display very high intramus- 
cular fat contents, fat mass and fat proportion, especially 
in males [22]. Previously, genetic modifier regions for the 
effect of the myostatin Mstn Cmpt ~ d!1Abc mutation on 
muscle mass have been identified on chromosomes 3, 5, 7, 
11, 16, and X in a cross between Comp9 and CAST/Ei 
lines [23]. 

Given that the BMMI lines were originally developed as 
supporting genetic models for livestock research, we were 
particularly interested in myostatin effects on intramuscu- 
lar fat content (IMF) and water holding capacity (WHC) 
on the genetic background of high muscularity. We also 
analysed the extent to which sex and the direction of the 
reciprocal cross impacted on the traits of interest. The 
latter could indicate parent-of-origin effects, where the 
impact on the phenotype can be different depending on 
the parent from which an allele was inherited. For ex- 
ample, the polar over-dominance caused by the ovine 
callipyge locus, where a hypermuscular phenotype only 
occurs if the mutated allele is inherited from the sire 
[24,25]. Parent-of-origin effects have also been described 
for body composition and fat-related traits in mice, pigs 
and cattie [26-29]. 



In addition to the relationship between muscle mass 
and meat quality traits, we were also interested in certain 
parameters of the muscle and whole body metabolism 
such as muscle glycogen and lactate contents, blood glu- 
cose levels, and the carcass pH-values. For this purpose, 
we present the correlations between these traits in the 
G3-population. The linkage study did not reveal genomic 
loci accounting for variation of those metabolic traits. 

Results and discussion 

Phenotypes 

As shown in Table 1, significant differences were found 
between the two parental lines. Averaged over both sexes, 
the Mstn Cmpt ~ dllAbc mutant BMMI866 animals showed 
42%, 42%, 99% and 94% higher values for body weight, 
lean mass, M. longissimus and M. quadriceps masses than 
the BMMI806, respectively. BMMI866 mice had 30% 
lower total fat percentage than BMMI806. The IMF 
contents of the M. longissimus and the M. quadriceps 
were 52% and 40% lower, respectively, as compared to the 
BMMI806 line. These data confirm the hypertrophic effect 
of the Mstn Cmpt ~ dllAbc mutation and its impact on fat 
accumulation [7]. Furthermore, fasting blood glucose 
levels of the BMMI866 line were 14% below the levels of 
BMMI806. The decreased glucose levels of BMMI866 
mice support the model of a metabolic shift towards the 
utilization of glucose as energy fuel if myostatin is not 
fully functional, as shown by experiments in cell cultures 
[30]. Regarding differences between the sexes, male 
BMMI866 mice had lower carcass pH-values after 1 hour 
and female BMMI866 mice showed lower carcass 
pH-values after 24 hours post-mortem compared to 
BMMI806 (p < 0.05). 

Comparing the Fi and G3-population with the parental 
lines, we found intermediate values for body weight, lean 
mass and muscle masses. Dominance of the BMMI806 
characteristics was indicated for total fat percentage in 
the Fi population. While the IMF values of both sexes of 
Fi animals and the males of the G 3 population were 
intermediate, in female G3 animals, the measurements of 
total fat percentage and IMF of the M. longissimus were 
similar to the lower values of the BMMI866 line. Fasting 
blood glucose levels indicated a dominance of the 
BMMI866 line alleles, whereas the BMMI806 phenotype 
was found prevalent for the higher carcass pH-values for 
1 hour post mortem. 

Correlations 

In the G3-population, high positive correlations in males 
and females were observed between muscle mass and total 
lean mass (r > 0.79) and negative correlations with body 
fat percentage (r < -0.46) and IMF (r < -0.53; Table 2). In 
females, muscle mass was also negatively correlated with 
WHC in M. quadriceps (r = -35), whilst in males, no 
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Table 1 Body traits of parental, F 1 and G 3 animals (means and standard deviations) 
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Trait 


Males 


Females 


Males 


Females 


Males 


Females 


Males 


Females 


Body weight (g) 


49.31 


40.35 


36.43 (4.69) 


26.77 (2.16) 


42.17 (344f°' bb 


34.38 


44.98 


36.69 


{Z.5b) 


U.4/J 








U- 1 z) 


[A 1 \aaa,bbb 

(4.2 I ) 


(4.63) 


Lean mass (g) 


a 1 on 
4 I .oy 


33. /y 


zy.3 i [z.Zy) 


Z3.0/ { i .yyj 


3Z.3U 


Zy.KjQ 


DO. 1 Z 


jU.O I 


(2.43)*** 


(2.12)*** 






(2 72) aoa 'kkk 


(2 03) aoa 'kk6 


(4.03) aaabbb 


(4.42) a<M ' bbb 


ML-Mass (g) 


1.33 (0.13)*** 


1.08 (0.1)*** 


0.67 (0.08) 


0.54 (0.05) 


0.89 (0.09) aaabbb 


0.75 (0.08) oaa ' bbb 


0.95 (022) oao ' bbb 


0.85 (023) aao ' bbb 


MQ-Mass (g) 


0.84 (0.05)*** 


0.67 (0.05)*** 


0.44 (0.04) 


0.34 (0.03) 


0.55 (0.05) aaabbb 


0.47 (0.03) oaa ' bbb 


0.61 (0.15r a ' bbb 


0.55 (0.15f°' bbb 


Fat mass (g) 


3.53 (1.67) 


3.50 (1.35) 


4.41 (1.46) 


3.13 (1.19) 


4.41 (0.98) a 


3.76 (0.94) 


5.00 (1.96)°° 


2.88 (1.39) a 


Fat mass (%) 


7.77 (3.65)*** 


9.36 (3.49) 


12.90 (3.52) 


11.41 (347) 


1 1 .98 (2.35) atM 


11.40 (2.31)° 


12.16 (4.50) aaa 


8.70 (4.3 1) b 


MF-ML (%) 


1.17 (0.23)*** 


1 .28 (0.23)*** 


2.64 (0.34) 


2.45 (0.49) 


1 .89 (0.32) aao ' bbb 


1 .50 (0.2f' bbb 


1 .69 ( 0 .39) oaa ' bbb 


1 .38 (0.33) bbb 


IMF-MQ (%) 


1.07 (0.15)*** 


1.11 (0.12)*** 


1.81 (0.30) 


1 .80 (0.29) 


1.46 (0.1 5f aa ' bbb 


1.38 (0.16f aa ' bbb 


1.34 (021) oaa ' bbb 


1.22 (0.18r' bbb 


WHC-ML (%) 


1 .00 (0.58) 


1 .03 (0.39) 


1 .09 (0.37) 


1.01 (0.39) 


0.90 (0.25) fa 


1 .08 (0.35) 


1 .24 (0.64)°° 


1.23 (0.41 f 


WHC-MQ (%) 


0.75 (0.21) 


0.69 (0.17) 


0.85 (0.33) 


0.81 (0.26) 


0.86 (0.1 If 


0.95 (0.33)°" 


0.93 (0.29)°° 


0.93 (0.23) aao ' b 


Glucose (mg/dL) 


105 (24)* 


88 (17)** 


120 (26) 


104 (17) 


108 (14) 


87 (12) bbb 


106 (13) bb 


90 (12) bbb 


Glycogen (mg/g) 


1.94 (0.75) 


2.24 (1.51) 


247 (0.76) 


2.45 (1 .38) 


1 .52 (0.52) b 


1 .73 (0.54) 


1 .37 (0.90) bb 


1.18 (0.76) b 


Lactate (mg/g) 


0.068 (0.008) 


0.091 (0.014) 


0.064 (0.003) 


0.070 (0.025) 


0.035 (0.010) ao ' bb 


0.035 (0.01 2) ab 


0.048 (0.017) oa ' bb 


0.045 (0.013) aa ' bb 


pH 1 hour 


6.67 (0.11)* 


6.63 (0.08) 


6.74 (0.09) 


6.71 (0.12) 






6.77 (0.1 5f a 


6.73 (0.15) aa 


pH 24 hours 


6.16 (0.10) 


6.15 (0.12)* 


6.12 (0.12) 


6.25 (0.12) 






6.12 (0.11) 


6.13 (0.09) bbb 



Trait abbreviations: IMF = intramuscular fat content, WHC = water holding capacity (measured as 'drip loss' in percent of muscle weight), ML = Musculus 
longissimus, MQ = Musculus quadriceps, Glucose = blood glucose level. Phenotypes of the parental lines are compared with each other and the Ft and G 3 , 
separated for sex. Asterisks indicate significant differences between the parental lines; letters indicate significant differences between the Ft or G 3 and BMMI866 
(") or BMMI806 (") lines. Number of asterisks and letters show levels of significance ("nwWM* p < 0 001 ( ,,/aa/bb p < Q m Mb p < Q Q5) 



correlation was found. Furthermore, in males, higher 
muscle mass was associated with lower pH-values 
24 hours post mortem (r < -0.16). Nevertheless, in 
both sexes negative correlations were evident between 
muscle glycogen contents and the pH-values of 1 and 
24 hours post mortem (r < -0.15), and between WHC and 
the pH-value after 24 hours (r < -0.14). The muscle glyco- 
gen content was associated with the body fat content 
(r < -0.17), but was not correlated with the blood glucose 
levels; higher muscle glycogen contents were found in lea- 
ner animals. 

The results regarding the carcass pH-value and 
WHC confirm studies in pigs, where a low pH-value 
negatively affects WHC and leads to high drip loss. 
The low WHC is likely the result of a developing acidic 
environment that affects the proteolysis of cell scaffold 
proteins [31,32]. In turn, the acidic environment is 
caused by high intramuscular glycogen contents that 
promote the development of high lactate levels via gly- 
colysis [31]. The negative correlation between muscle 
glycogen and body fat content as well as IMF of the 
M. longissimus in males supports the model of a meta- 
bolic shift towards glycolysis. The non-significant correl- 
ation with fasting blood glucose levels is likely due to 
the short fasting period of two hours before measure- 
ments leading to a high variability. 



QTL analysis 

Due to the relative small number of animals in our 
G3-population, we could only detect regions with larger 
effects. Two genome-wide significant QTL were identi- 
fied (Figure 1). The first one was located on chromosome 
(Chr) 1 at 54 Mb (47-59 Mb) in a region that contains the 
previously described myostatin Mstn Cmpt ' dllAhc mutation, 
located at 53 Mb (Figure 2). In the G 3 -population, the 
BMMI866 allele at this locus caused additive effects of 
1.25 g for body weight and 1.35 g for lean mass. It also 
increased the masses of the M. longissimus and M. quadri- 
ceps muscles by 0.07 g and 0.04 g, respectively (Table 3). 
The body fat percentage was decreased by the BMMI866 
allele by 1.15 points. A similar effect was observed for 
IMF where the BMMI866 allele had a decreasing effect of 
0.13 and 0.07 percentage points for IMF values of the 
M. longissimus and M. quadriceps, respectively (Figure 3a, 
Table 3). 

An additional QTL was identified on Chrl6 at 94 Mb 
(86-94 Mb; Figure 2) that decreased the total fat mass by 
0.55 g in heterozygous animals (Figure 3a, Table 3). This is 
a negative heterosis effect where heterozygous animals have 
significantiy lower t values than either homozygous class. 
This significant negative dominance for fat mass at the Chr. 
16 QTL is interesting because it is in contrast to the pheno- 
typic observation that male Fl and G3 individuals tended 



Table 2 Spearman's correlation coefficients between different traits in the G 3 population 
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Weight (g) 


Mass (g) 


(g) 


(g) 
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(%) 
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(mg/dL) 


(%) 


(%) 


(mg/g) 


(mg/g) 


1 hour 


24 hour: 


Body weight (g) 
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U.00 
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— U. 1 z . 


—U.J 


-U.Zj 


u.u 


—u.u I 
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n 1 
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Lean mass (g) 


u.yj 




U.oU 


n 7G*** 

u./y 


— U.zV 


— u.du 


— U.3U 


— U.4y 


— U.Uo 


U.U I 


— U.U4 


U. 1 5 


—U.Uo 


—U.U/ 


— U. I / 


ML-Mass (g) 


0.82*** 


0.87*** 




0.91*** 


-0.32*** 


-0.51*** 


-0.65*** 


-0.69*** 


-0.01 


0.04 


-0.04 


0.20** 


-0.02 


-0.12 


-0.19** 


MQ-Mass (g) 


0.82*** 


0.87*** 


0.96*** 




-0.28*** 


-0.46*** 


-0.66*** 


-0.64*** 


-0.03 


0.05 


-0.1 


0.11 


-0.04 


-0.12. 


-0.16* 


Fat mass (g) 


-0.10 


-0.32*** 


-0.28** 


-0.29** 




0.96*** 


0.52*** 


0.62*** 


0.03 


-0.01 


-0.03 


-0.13 . 


0.09 


-0.14* 


-0.11 


Fat mass {%) 


-0.34*** 


-0.54*** 


-0.48*** 


-0.49*** 


0.96*** 




0.60*** 


0.70*** 


0.04 


-0.02 


-0.02 


-0.17* 


0.08 


-0.1 


-0.05 


IMF-ML (%) 


-0.51*** 


-0.60*** 


-0.64*** 


-0.66*** 


0.45*** 


0.56*** 




0.67*** 


0.05 


0.05 


0.03 


-0.15* 


-0.01 


0.03 


0.08 


IMF-MQ (g) 


-0.35*** 


-0.45*** 


-0.57*** 


-0.53*** 


0.37*** 


0.46*** 


0.52*** 




0.03 


0.03 


0.09 


-0.12 


-0.09 


0 


0.08 


Glucose (mg/dL) 


0.05 


0.05 


-0.08 


-0.01 


0.01 


0.01 


0.13 


0.11 




-0.01 


0.01 


0 


-0.03 


-0.05 


-0.05 


WHC-ML (%) 


-0.16 . 


-0.13 


-0.15 


-0.09 


-0.16 . 


-0.09 


-0.05 


0.03 


0.15 




0.35*** 


0.15* 


-0.02 


0.08 


-0.21** 


WHC-MQ (%) 


-0.38*** 


-0.39*** 


-0.35*** 


-0.35*** 


0.04 


0.13 


0.09 


0.1 


0.04 


0.42*** 




0.10 


0.10 


0.11 


-0.14 . 


Glycogen (mg/g) 


0.09 


0.10 


0.13 


0.12 


-0.24* 


-0.23* 


-0.04 


-0.02 


0.14 


0.05 


0.03 




-0.05 


-0.16* 


-0.27*** 


Lactate (mg/g) 


0.04 


0.06 


0.10 


0.08 


-0.05 


-0.07 


-0.14 


-0.20* 


0.06 


0.14 


0.15 


0.07 




-0.17* 


-0.15* 


pH, 1 hours 


0.05 


0.07 


-0.11 


-0.12 


0.01 


-0.01 


0.12 


0.08 


0.01 


-0.02 


-0.06 


-0.20* 


-0.17 . 




0.18** 


pH, 24 hours 


0.03 


0 


-0.08 


-0.07 


-0.01 


-0.01 


0.04 


0.01 


0.01 


-0.23* 


-0.16 . 


-0.31*** 


-0.23* 


0.23* 




Females 

































Trait abbreviations: IMF = intramuscular fat content, WHC = water holding capacity, measured as 'drip loss' in percent of muscle weight, ML = M. longissimus, MQ = M. quadriceps. Glucose = blood glucose level; Asterisks 
indicate the significance of correlations: *** p < 0.001 , ** p < 0.01 , * p< 0.05, . p < 0.1 0. 



Karst et al. BMC Genomics 2013, 14:16 
http://www.biomedcentral.com/1471-2164/14/16 



Page 5 of 1 2 



0,0 P i centromere 



1 



rs31991963 || |^ |g || || |g J? 



centromere 0.0- 



rs36990S9 
rs1 3476603 
■rs274 10285 
rs27369144 



rs27249426 
■rs13476877 1 



centromere 0,0 " 



rs3036375S 88.6- 
95,0- 

■rs30709106 

117,3- 



rs27817407 30,7- 



154.6 -JL" 



rs3696703 152.5- 



•rs134 78166 
■ rs33E 15624 



- rs3719351 
-rs29819195 



■ rs32488861 jg 3 . 
-rs32154092 37' . 



-rs31 697976 

-rs31 045562 
•rs31823507 , 
- rs32325529 1 



- centromere 0.0 — 

-rs3683183 
-rs1 34801 21 
-rs3691977 



■ rs3696502 

■ rs33507556 



■rs32760317 
rs13462977 124,1 



- rs29359796 

- rs2937865C 



Chr12 

— centromere 



-rs29391262 
- rs27059666 



-cenlromere 0,0 
-re3719002 8,9 

- rs13481720 

24,6 

-rs13481770 38,1 

- rs6249046 4g ? 

- rs8243445 

59 5 

-rs3703511 

73,4 



- rs13481961 
-rs3719361 , 



-rs30953521 8 6- 

-rs3!333715 2 5 6- 

- rs30916926 34 4 " 
-rs31327249 45 3 ~ 

- rs30252928 
-rs30742840 76 3 . 



■ rs32256782 2 
rs31656820 
rs31996797 4 



■rs32460616 93,6- 
■rs32007840 98 3- 



-rs4167212 
-rs4179051 



-rs4197590 
- rs420S997 



- rs33583385 63,0 

- ts33082964 
-rs33717365 q 0ft 



— centromere 

- rs3700209 



- rs30563254 
TS30983093 



Figure 1 MapChart plot of 138 reference single nucleotide polymorphisms used in this study. Positions are given in Mb (NCBI Build 37) 
Bars indicate LOD support intervals of identified significant QTL 



to have a higher fat mass than the mean of to the parental 
lines (Table 1). The general positive heterosis of the off- 
spring lines points to other loci with effects on fat mass that 
could not be detected in this study. The QTL on chromo- 
some 16 showed no linkage with IMF, which indicates a 
specific role in adipose tissue for this region in mice. An 
interaction between the Chrl and the Chrl6 locus was not 
observed, and thus our data provides evidence that the re- 
gion on Chrl6 is an independent regulator of total fat mass. 
Up to now only one marker in that region, the D16Mit51 
at 93 Mb, was proposed provisionally by Srivastava et al. 
(2006) as causative for variation in total fat mass in mice. 

The QTL for fat mass on Chrl6 is located within a 
region from 75-98 Mb that is syntenic with a region 
in pigs on Sus scrofa chromosome (SSC) 13 (129-145 Mb). 
According to the animal genome database (animalgenome. 
org) this region contains QTL for marbling, back fat and 
for IMF [32-35]. Dominance effects for fat mass, as they 
were found in our study, were not significant for the syn- 
tenic region in pigs, however a positive dominance was 
found for ham weight [36-38]. Thus, this syntenic locus 
could possibly have a negative dominance effect on fat 
mass because muscle weight and fat weight are often 
negatively correlated. Furthermore, similar effects were 
observed for the syntenic region in cattle. In cattle, this 
region is divided into two parts which reside on chromo- 
some 1 (BTA1) between 0-23 Mb and 141-155 Mb. 
Within those two regions on BTA1, fat related QTL were 



mapped for milk fat yield, marbling and fat thickness at the 
12th rib [39-42]. A positive but non-significant dominance 
effect on adjusted subcutaneous fat thickness was reported 
in cattle [40]. Finally, the syntenic region in the human 
genome resides on chromosome 21 at 15-43 Mb. The 
Human Obesity Gene Map counts four obesity-associated 
loci in this syntenic region [43,44]. Assuming that the same 
gene is underlying the variation across species, further 
refinement of the murine Chrl6 regions seems promising 
and our study adds to the confirmation of mice as a model 
animal to study meat quality traits in livestock or muscular 
disorders in humans. 

As we did not find any other loci affecting muscle mass 
apart from the myostatin locus, we assume that several or 
many small effect alleles that contribute to these traits and 
that have been shown in other crosses [45] are hidden in 
our study behind the strong effect of myostatin, in par- 
ticular in the hypermuscular line BMMI806. Another 
reason could be that the statistical power was not suffi- 
cient in the examined population. 

Effects of sex and direction-of-cross 

As shown in Table 3, sex affected all traits for all identi- 
fied QTL. Males significantly increased body weight, 
lean mass, and the masses of the M. longissimus and 
M. quadriceps by 8.43 g, 5.74 g, 0.10 g and 0.07 g, res- 
pectively. In males, the total fat mass and the total fat 
percentage were also raised by 1.98 g and 2.97 points, 
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respectively. Sex also affected IMF values of the two 
examined muscles, which were increased by 0.28 points 
and 0.41 points for the M. longissimus and the M. quad- 
riceps in males, respectively. The increase in the lean 
mass related traits was expected as a sex difference; 
however, the increase in fat mass percentage in male 
mice was a new observation. An explanation could be 
that male BMMI806 mice display very high fat values 
(Table 1) and the high-fat phenotype in G 3 -males 
is thought to be inherited by the BMMI806 alleles of 
both QTLs. 

Interestingly, for the myostatin locus, we also found 
significant effects of the direction-of-cross and an inter- 
action between sex and the direction-of-cross on IMF of 
the M. quadriceps (LOD = 13.62). Animals derived from 
mating a female BMMI806 with a male BMMI866 showed 
higher body weight, lean mass and muscle mass at the 
Chrl locus in comparison to animals of the reciprocal 
cross. While total fat mass, fat percentage and IMF of the 
M. longissimus were not significantiy affected, we observed 
a large increase of IMF in the M. quadriceps by 0.29 
points in male G 3 -animals that were derived from a 
mating between a female BMMI866 and a male BMMI806 
(Figure 3b,c). This means that in male G 3 -animals 



descending from a BMMI806 grand-grandfather, myosta- 
tin expressed less of its IMF decreasing effect. Further- 
more, the effects of myostatin were generally larger in 
G3-animals descending from a male BMMI866 parent of 
the initial cross. 

As a likely result of the recent selection and inbreed- 
ing process, SNP data of the mitochondrial DNA and 
the Y-chromosome did not differ between the two par- 
ental lines. Therefore, we assumed that mitochondria 
and Y-chromosomes did not account for trait differences 
between reciprocal crosses. Nevertheless, we cannot 
completely exclude hidden mutations derived during the 
breeding history. Such hidden mutations do not play a 
role for mitochondria, but might affect Y-chromosomal 
effects. However, so far, no Y-chromosomal effects on 
traits that we analysed in this study have been identified. 
Therefore, the phenotypic differences between the recip- 
rocal crosses are likely parent-of-origin effects. 

The high IMF values for the QTL on Chrl, for G 3 -males 
of the cross between a BMMI866 female and a BMMI806 
male, were very similar to males of the parental line 
BMMI806. This was not observed in the reciprocal cross 
and therefore, might result from parent-of-origin effects. 
Since evidence is growing that parent-of-origin effects 



Karst et al. BMC Genomics 2013, 14:16 
http://www.biomedcentral.com/1471-2164/14/16 



Page 7 of 12 



Table 3 QTL, sex and direction-of-cross effects (p < 0.001) in the G 3 population at age of 10 weeks 



Trait 


QTL/Effect 1 


Chr 


Mb 2 


CI 3 


Marker 4 


LOD 5 


a (s.e.) 6 


d (s.e.) 6 


% G 3 var 7 


Body weight (g) 


sex 










83.01 


8.43 (0.31) 




68.35 




pgm 










6.44 


1.67 (0.3) 




2.81 




1<a>27.0 


1 


54 


47-59 


rs31991963 


8.93 


-1.25 (0.21) 


-1.03 (0.31) 


3.97 


Lean mass (g) 


sex 










62.43 


5.74 (0.27) 




55.22 




pgm 










7.78 


1 .59 (0.26) 




4.42 




1<a>27.0 


1 


54 


47-55 


rs31991963 


13.91 


-1.35 (0.18) 


-1.2 (0.26) 


8.27 


ML-Mass (g) 


sex 










14.30 


0.1 (0.01) 




15.31 




pgm 










3.43 


0.04 (0.01) 




3.38 




1(327.0 


1 


54 


47-59 


rs31991963 


17.26 


-0.07 (0.01) 


-0.07 (0.01) 


18.90 


MQ-Mass (g) 


sex 










8.68 


0.07 (0.01) 




10.86 




pgm 










2.50 


0.03 (0.01) 




2.98 




1(327.0 


1 


54 


47-55 


rs31991963 


7.17 


-0.04 (0.01) 


-0.04 (0.01) 


8.88 


Fat mass (g) 


sex 










35.31 


1.98 (0.15) 




40.23 




16(355.0 


16 


94 


86-94 


rs3143713 


3.57 


-0.07 (0.12) 


-0.55 (0.15) 


3.16 


Fat mass {%) 


sex 










15.49 


2.97 (0.36) 




18.85 




1(327.0 




54 


47-59 


rs31991963 


6.32 


1 .03 (0.24) 


1.15 (0.35) 


7.17 


MF-ML (%) 


sex 










20.04 


0.28 (0.03) 




22.55 




1(327.0 


1 


54 


47-55 


rs31991963 


10.15 


0.13 (0.02) 


0.06 (0.03) 


10.56 


IMF-MQ (%) 


sex 










50.46 


0.41 (0.02) 




37.56 




pgm 










39.19 


-0.01 (0.03) 




26.50 




1 (326.0 


1 


51 


47-59 


rs3 199 1963 


10.05 


0.07 (0.01) 


0.00 (0.02) 


5.36 




sex:pgm 










16.08 


-0.29 (0.03) 




9.00 



1 QTL positions are given as chromosome number, @, and the genetic distance from the centromer in cM. 

2 Most likely chromosomal location given as Mb-position. 

3 1-LOD support interval in Mb. 

4 Marker closest to the chromosomal position with the highest LOD score. 

5 LOD scores from the full model estimations. 

6 additive (a) and dominance (d) effect and their standard errors (s.e.) determined with the non-transformed raw trait values, therefore given in the respective 
unit; the direction of a and d, respectively given as the effect of the BMMI866-allele. 

7 G 3 phenotypic variance (%) explained by the QTL; QTL effect given as the reduction of the residual sum of squares fitting 1 vs. 0 QTL. 



could exist over several generations [46] our findings 
suggest epigenetic patterns, which alter the gene expres- 
sion of the chromosomes according to the ancestor they 
were inherited from. Whether the myostatin locus alone 
underlies genetic imprinting, or additional loci that did 
not directly affect the analysed traits, still has to be tested 
in a genome-wide screen. 

Recently, we have identified bidirectional parent-of-origin 
effects for muscle glycogen content and glycolytic potential 
as well as body weight in a population between two other 
hypermuscular mouse lines BFMI806 x BMMI816 [47]. 
The imprinted loci were not discovered by a genome-wide 
QTL search directly for additive and dominance effects. 

In sheep with the callipyge mutation, the callipyge 
phenotype was not expressed in homozygous mutated 
animals [24,25]. It was suggested that the maternal 
inactivation of the locus also negatively affected the 
paternal expression, but a paternally inherited mutated 
locus escapes the inactivation in the presence of a 



maternally inherited wild-type allele and leads to the cal- 
lipygous phenotype. Phenomenon like this have their 
molecular origin in a DNA methylation 'reset' that occurs 
during gametogenesis, the fusion of oocyte and spermato- 
cyte in the zygote and the following early embryogenesis 
[48,49]. As seen in the callipyge sheep example, small 
genetic differences can affect the reprogramming pro- 
cesses and lead to different gene expression patterns and 
phenotypes. 

Conclusions 

The strong effects of the myostatin mutation 
Mstn Cmpt ' dllAhc on high muscle and low fat mass were 
confirmed, even in a background of high muscularity in 
a cross between two hypermuscular mouse lines. The 
data provides evidence for sex and very likely parent-of- 
origin effects modifying the direct effects of the myosta- 
tin locus on muscle mass and IMF in M. quadriceps. 
However, further tests are necessary to confirm imprinting 
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(See figure on previous page.) 

Figure 3 Effect plots showing means and standard errors of the three genotype-classes at significant QTL positions, a) QTL significant 
under 95% and b) sex and direction-of-cross effects; c) directions-of-cross effects for the three genotypes for lean mass and IMF of the M. 
quadriceps; 866 = homozygous BMMI866-allele, 806 = homozygous BMMI806-allele, H = heterozygous animals; 806f x 866m = G 3 -animals 
originating from mating a female BMMI806 and a male BMMI866, 866f x 806m G 3 -animals originating from mating a female BMMI866 and a male 
BMMI806 in the parental generations. 



effects. In addition to the myostatin locus, we have shown 
that a region on Chr 16 affecting fat mass is conserved 
between mice, pig and cattle showing effects on body fat, 
marbling and milk fat. Furthermore, many reported obser- 
vations and correlations of the post mortem physiology of 
livestock were reproduced in our mice. 

The results obtained in this study help to justify the 
mouse model to study economically relevant livestock 
traits. New information from the mouse model helps to 
describe variation in meat quantity and quality and thus 
supports research in livestock. 



Methods 

Animals 

The Berlin Muscle Mouse population had been long 
term selected for high body weight and high muscle 
mass to reflect the selective mechanisms in livestock 
breeding. Founder animals of the Berlin Muscle Mouse 
(BMM) population were originally purchased in several 
pet shops in Berlin, Germany. The selection process 
comprised several distinct phases. The beginning of the 
selection process constituted a phase of 23 generations 
of selection for high protein content of the carcass at the 
age of 60 days. Protein content was determined by 
chemical analyses. In a second phase, mice were selected 
for high body weight and low fat content at 42 days for 
10 generations. Afterwards, mice were monitored for 
high muscularity by palpating and mice with highest 
muscularity on a scale between 1 and 5 were selected for 
the next generation. As a result of selection and likely 
the occurrence of natural mutation(s) over 25 genera- 
tions, a mouse population with a high muscular pheno- 
type had been generated. A high compact sub-line was 
perpetuated through random mating of selected animals 
[21]. Sequencing of the myostatin gene in this line 
revealed a 12 bp deletion [20] leading to a loss of func- 
tion. Scale and weight based selection continued for 
another 28 generations. After 86 generations of selec- 
tion, full-sibs with distinct phenotypes were mated. 
These founder animals became the basis of seven Berlin 
Muscle Mouse inbred lines (BMMI). Four lines display 
the Mstn Cmpt ' dllAbc mutation and three are wild type. In 
this study we used the lines BMMI866, which carries the 
described Mstn Cmpt ' dllAbc mutation, and BMMI806, 
which is hypermuscular but wild type for myostatin and 
shows high IMF. At the time of setting up the crossbred 



experiment, the lines were in generation 20 (BMMI866) 
and 21 (BMMI806) of inbreeding. 

Pedigree structure 

Two pairs of full sibs of the Berlin Muscle Mouse inbred 
lines BMMI866 and BMMI806 were crossed reciprocally 
to generate Fl, F2 and G3-intercross populations. F 2 -ani- 
mals within every direction-of-cross were randomly 
mated under avoidance of sibling-mating [50], to gener- 
ate 308 G3-animals. Altogether, there were 195 males 
and 113 females. From these 308 animals, 97 males and 
53 females came from the cross between a female 
BMMI866 and a male BMMI806 mouse, while the recip- 
rocal cross between a male BMMI866 and a female 
BMMI806 mouse consisted of 98 males and 60 females. 
Using this design, we would expect different mitochon- 
dria and Y-chromosomes segregation in the reciprocal 
crosses. Since high-density SNP data of the inbred 
strains BMMI866 and BMMI806 provided evidence for 
no variation in both mitochondria and Y-chromosomes, 
we assume that the two parental mouse strains have the 
same mitochondria and Y-chromosomes as a result of 
their breeding history (See also under Genotyping). 

Husbandry and feeding conditions 

The mice were treated in accordance to, and all experimen- 
tal protocols were approved by, the German Animal 
Welfare Authorities (approval no. G0405/08). The animals 
were maintained under conventional conditions at 22 ± 2°C 
and controlled lightning with a 12:12 hour lightdark cycle. 
They were kept in groups of two to four animals of the 
same sex per Macrolon cage and given ad libitum access to 
food and water. Until the age of 70 days, the animals were 
fed a standard breeding diet (Altromin standard breeding 
diet no. 1314 TPF| Lage, Germany). This diet contained 
27.0% crude protein, 5.0% crude fat, 4.5% crude fibre, 6.5% 
crude ash, 50.5% nitrogen free extract (starch and sugar), 
vitamins, trace elements and minerals (2988 kcal/kg me- 
tabolizable energy; thereof 27.0% energy from proteins, 
13.0% from fat and 60.0% from carbohydrates). 

Phenotypic measures 

After a fasting period of two hours, 71 day-old mice 
were anaesthetised under isoflurane and decapitated. 
The Musculus longissimus (ML) and Musculus quadri- 
ceps (MQ) were dissected and weighed. The summed 
muscle weights of the left and right M. longissimus and 
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the left and right M. quadriceps were recorded as muscle 
mass (MM). The right muscles were immediately frozen 
in liquid nitrogen and subsequently stored at -80°C. The 
left muscles were cooled down to 6°C for one hour, subse- 
quently frozen at -20°C and stored at the same tem- 
perature until WHC and IMF were measured. Carcasses 
were stored at 6°C and pH-values were taken within the 
M. biceps femoris at one and 24 hours post mortem (ebro 
PHT 810, Ingolstadt, Germany). For the determination of 
WHC, frozen muscles were thawed and stored at 6°C for 
24 hours. Muscles were then centrifuged for 60 sec at 
604 x g in invitek® 1.5 ml receiver tubes with filter inlays 
to collect the tissue fluid that was not held from the muscle 
(Eppendorf Minispin, Hamburg, Germany; Invitek, Berlin, 
Germany). The ratio of lost tissue fluid to tissue mass was 
designated as 'drip loss'. IMF was measured as percen- 
tage of muscle weight using nuclear magnetic resonance 
technology (SMART Trac System, CEM, Kamp-Lintfort, 
Germany) [51]. Body weights were recorded weekly. At ten 
weeks, total fat and lean masses were measured by quanti- 
tative magnetic resonance (QMR) analysis, using the 
EchoMRI whole body composition analyser (Echo Medical 
Systems, Houston, Texas, USA) [22,52]. Body fat percent- 
age is the percentage of total fat mass calculated from the 
sum of total fat and lean mass. Blood glucose levels were 
measured after two hours of fasting, prior to dissection at 
ten weeks (Bayer 'Contour! Leverkusen, Germany). Muscle 
glycogen content was determined colorimetrically in the 
right M. longissimus (GOD/PAP method 'Glucose liquico- 
lor' by Human, Wiesbaden, Germany) as suggested by 
Barham and Trinder (1972). Lactate contents were deter- 
mined colorimetrically in the right M. longissimus using 
the Lactate Assay Kit by Techung Lee [53]. 

Genotyping 

Parental BMMI lines were genotyped with the Mouse- 
Diversity-Array [54] comprising 623,124 single-nucleotide 
polymorphisms (SNPs). The SNP information provided 
evidence for allele fixation of 98.1% and 98.3% in lines 
BMMI866 and BMMI806 after 20 and 21 generations of 
inbreeding, respectively. Both lines differed from each 
other by 4.8% at the SNP level. SNPs of the mitochondria 
and the Y-chromosome were not different between the 
lines. Using the information on diverse genomic regions 
between BMMI806 and BMMI866, 138 informative SNP 
markers covering all chromosomes (except Y) in an aver- 
age distance of 24.8 Mb were selected for genotyping the 
parents, F 2 and the G 3 -animals (Figure 1). Regions larger 
than 10 Mb that did not differ in SNP-alleles between 
parental lines could not be included into the linkage ana- 
lysis. We assume that these regions have the same ances- 
tral origin as a result of the breeding history. Genotyping 
was done at KBiosciences (Hoddesdon, U.K.). The phys- 
ical map was converted into the genetic map using the 



'Mouse Map Converter' software from The Jackson 
Laboratory [55]. 

QTL and statistical analyses 

For the QTL analysis, 308 G3-animals were used. Analyses 
for single QTL detection and detection of interacting 
QTL were performed using R/qtl [56]. Unlike F 2 , the 
G3-animals may be unequally related to each other. Ignor- 
ing the unequal relatedness may result in a serious infla- 
tion of false positive rates. Therefore, we calculated 
environmental residuals for both sexes and each of the 
two crosses and corrected for litter size and genomic kin- 
ship. This part of the analysis was carried out with a 
Genome-wide Rapid Analysis using Mixed Models And 
Mixed Models and Regression (GRAMMAR) as imple- 
mented in GenABLE [57-60]. The environmental residuals 
were used to generate corrected phenotypes, which were 
applied in the actual QTL analysis. Direction-of-cross and 
sex were included as additive covariates in the model and 
used as interactive covariates to test their effects on QTL. 
All phenotypes were log-transformed to obtain normal 
distribution. Trait-specific significance thresholds were 
estimated with 1000 permutations [61]. On average, the 
thresholds for genome-wide significance and suggestive 
linkage results were LOD = 3.59 (p < 0.05) and LOD = 3.19 
(p < 0.10), respectively. Genome-wide significant QTL 
(p < 0.05) were included into a mixed model for each trait 
to calculate the respective trait variance in the G 3 -popula- 
tion. Factors being significant under p < 0.001 were kept in 
the model. The genotype of the myostatin locus was 
accounted for by including it as a covariate in scans for 
additional QTL. Basic statistics were performed using the 
SAS software package (SAS System for Windows, Release 
9.2). QTL support intervals were determined by bootstrap- 
ping as implemented in R/qtl [60]. 
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